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Abstract 

We investigate the decays D° 'ir~i'^u and D° K~£'^u, where ^ is e or /j,, us- 
ing approximately 7 fb~^ of data collected with the CLEO III detector. We find Rq = 
B{D^ ■K-e'^v)/B{D^ R-e+u) = 0.082 ± 0.006 ± 0.005. Fits to the kinematic distributions of 
the data provide parameters describing the form factor of each mode. Combining the form factor 
results and i?o gives |/^(0)|2|F,,|V|/f (0)|2|K,|2 = 0.038lO:°K°°5. 

PACS numbers: 13.20.Fc, 14.40.Lb, 12.38.Qk 
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The quark mixing parameters are fundamental constants of the weak interaction. Mea- 
suring them also tests the unitarity of the quark mixing (CKM) matrix, which is sensitive 
to as yet undiscovered particles and interactions. Semileptonic decays have provided most 
quark coupling data. For these decays, the strong interaction binding effects, parameterized 
by form factors, are simplest to calculate; nonetheless, even here, form factor uncertainties 
can dominate the experimental uncertainties 0. 

We present a study of the decays — > 7r~£"'"z/ and D° K'i^u, where i = e or /i. 
Charge conjugate modes are implied throughout this paper. We measure the ratio of their 
branching fractions, Rq = B{D^ 'K^e^v)/B{D^ K^e^u), and, for the first time for 
-K^l^v, parameters describing their form factors. The study of the 7r~£+z/ form 

factor is particularly interesting because it tests predictions for that of the closely related 
decay n~i~^h', which provides \Vub\- 

In the limit (m^/mc)^ = 0, where me and rric are the lepton and charm quark masses, the 
differential partial widths for 7r~£+z/ and — > K~l^v, in terms of the form factor 

are 

Here h is n or K, and is the invariant mass squared of the lepton- neutrino system and 
ranges from rnj to 2.98(1.88) GeV^ for TT~{K~)i~^i>. To reduce the form factor sensi- 

tivity of -Ro and determine the distributions, the yields are extracted in bins of g^. 

We use e~^e~ cc events collected at and just below the T(45') resonance with the 
CLEO III detector 0. We use only runs with good lepton identification, which leads to 
slightly different, but overlapping, datasets for the electron and muon modes with integrated 
luminosities of 6.7 and 8.0 fb~^ respectively. 

A major challenge for this analysis is the contamination of the D° ti~1^v sample by 
K~i'^i' decays, which are about a factor of 10 more common. The use of a Ring- 
Imaging Cherenkov detector (RICH) and specific ionization in the drift chamber {dE/dx) 
reduces this contamination dramatically by distinguishing K from vr mesons. The resulting 
efficiency and misidentification probability suppress misidentified —>■ K^i^u decays to 
15% of the D° 7i-£+u signal. 

The analysis also benefits from the hermeticity of the detector, which enables us to 
substitute the missing momentum vector of each event for the neutrino momentum. Within 
the active region, which covers 93% of the solid angle, we accept photons with energies above 
50 MeV and detect over 92% of charged particles with momentum above 75 MeV. 

candidates are reconstructed from lepton, hadron (vr or K), and neutrino combi- 
nations. Electron candidates have pe > 0.6 GeV, lie within the barrel of the detector 
(I cos6' I < 0.8, where 9 is the angle between the track and the beam) and have the expected 
calorimeter, RICH, and dE/dx signals. Muon candidates have > 1.5 GeV, lie within 
0.1 < I cos 6' I < 0.6, penetrate at least 5 interaction lengths of material, and have the ex- 
pected energy deposit in the calorimeter. The hadron (h) must have electric charge opposite 
that of the lepton and satisfy strong K or n identification requirements. The missing mo- 
mentum of the event (pmiss) provides the first estimate of the neutrino momentum: it is the 
negative of the net momentum of all charged particles and calorimeter showers (treated as 
photons) that are not associated with a track. We require rrihiu > 1-6 GeV. 

To improve the neutrino momentum resolution, we impose the mass constraint rrihiu = 
rriDO (we use PDG Q masses throughout). An ellipsoid of neutrino momenta satisfies this 
requirement; we take the momentum that lies in the plane defined by Pmiss and phi (hadron- 
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lepton momentum), and has the smallest vector difference from the missing momentum. This 
procedure reduces the full width at half-maximum of the neutrino momentum resolution from 
0.8 to 0.45 GeV. 

The kaons and pions from properly reconstructed decays tend to have higher momentum 
than those in backgrounds, so we demand ph > 0.5 GeV. We further require 0.6 < rrihi < 1.85 
GeV and \phiu\ > 2.25 GeV (low and middle bins) and \phiu\ > 3.0 GeV (high bin). 
Semileptonic B decays and most Bhabha, 2-photon, and t~^t~ backgrounds are suppressed 
by imposing 0.20 < R2 < 0.85, where R2 is the ratio of Fox- Wolfram moments j^. Bhabha 
and 2-photon events are also suppressed by demanding | cos ^thrust | < 0.8 for events in which 
a candidate e"*" (e~) lies in the hemisphere opposite the incident e+ (e~) beam. Here 6'thrust 
is the angle between the thrust axis of the event and the beam. 

We require that all candidates come from the decay D*~^ D^tt^. We reconstruct 
the D*^ by pairing a pion with the appropriate charge (the "soft" pion, tTs) with the 
candidate and then computing the mass difference between the D*~^ and the candidates. 
Am = rrihiuT^^ — rrihii,- The signal peaks in the region Am < 0.16 GeV (the "signal region") 
with a root-mean-square width of about 10 MeV. We use the Am distribution to extract 
the yields. 

About half the background in the signal region is composed of candidates in which the 
7r+ comes from a D*~^ decay. This background is troublesome because it peaks in Am, 
albeit more broadly than the signal, and, for ti~C.^v, is comparable in size to the 

signal. A Monte Carlo simulation [5J shows that most of this "peaking" background in 
the — s> TT^f'^z/ sample comes from K~t^v decays in which the K is mistaken for 

a vr (8%), or from candidates in which a lepton from K~C.^v (44%), K*~l^u 

(32%), —* p^i^u (9%), or nonresonant -D° {Ktt)~1^v (2%) is paired with a random 
pion or one from the same decay. The remaining half of the background does not peak 
because the tTs is not from a D* decay (the "false- tTs" background). For the more common 
K~l^v mode, the ratios of both the peaking and false-vTs background to signal are 
smaller by a factor of ten. The peaking background comes primarily from K*~i~^i' 
(66%), (Kn)-i+u nonresonant (6%), and D° n-i+u (4%). 

We divide the data into three bins: [0, 0.75] (bin 1), [0.75, 1.5] (bin 2), and > 1.5 
GeV^ (bin 3). The bin size is guided by our resolution of 0.4 GeV^. To calculate 
for —>■ K^i^u, we use in place of mx so that the K^i^u yield in each bin 

corresponds to the —>■ K~i~^i' background in the same n^i^u bin. 

The yield in each bin for each of the modes, K'e^u, K^fi^u, 

-K^e^v, and — >■ 7r~n~^i>, is determined from a fit to the Am distribution. The 
Monte Carlo simulation provides the Am distributions of the signal and backgrounds. 
The D° K~t^v samples are fit first. The two free parameters in these fits are the nor- 
malizations of the K~£~^i' simulated signal and of the false-TTg background relative to 
the data. Since the fit can only weakly distinguish the signal from the peaking backgrounds, 
we fix their ratio to the value predicted by the Monte Carlo simulation. (This assumption is 
investigated in the section on systematic uncertainties.) Then the D° ti^I^v samples are 
fit. The normalization of — > K~i^v from the — >■ K~i~^i> fits sets the normalization of 
the peaking background in the — > 7r~£~^i> fits. The two free parameters in these fits are 
the normalizations of the rt'l^u signal and of the false-TTg background. The electron 
mode fits and their confidence levels are shown in Figure ^ The muon fits are similar, but 
with smaller sample sizes because of the muon momentum and angular restrictions. To test 
for sensitivity to the details of the fitting shape, we reanalyze the —>■ K~e^v sample us- 
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FIG. 1: The fits to the Am distributions for — > K'e'^v (left) and — > TT~e'^i' (right) and their 
confidence levels (C.L.). The data (points) are superimposed on the sum of the normalized sim- 
ulated signal (peaked histogram), peaking background (dark histogram) and false- tTs background 
(broad histogram). 



ing Am calculated with the uncorrected neutrino momentum and also using Am calculated 
omitting the neutrino. We see no significant variation in the results. 

An efficiency matrix relates the number of decays produced in each bin (the efficiency 
corrected yields) to the number detected in each bin. Calculated using a Monte Carlo 
simulation, it accounts for both reconstruction efficiency and event migration across bins. 
The average reconstruction efficiency for vr^e+z/, not including the D*^ D^'^t 

branching fraction, is about 11%, and 30% of reconstructed events migrate from their true 

bin into another bin. For — >■ K~e~^i', the migration is somewhat greater because 
we use the pion mass to compute q^. Efficiencies are lower by a factor of four (six) for 

We sum the efficiency corrected yields over bins to find Rq^ = 0.085 ± 0.006 ± 0.006 
and Rq^ = 0.074 ± 0.012 ± 0.006 for the electron and muon modes, respectively, where the 
first uncertainty is statistical, and the second is systematic and is described below. We 
then compute the normalized distributions, which are defined as the fraction of the total 
corrected yield in each g^ bin (since the D* production rate is undetermined). They are 
shown in Figure El and Table IH The results combine the electron and muon modes after 
correcting the muon modes for their reduced phase space. The correlations between the g^ 
bins are pfg = -0.81, pfg = 0.18, p| = -0.72 and p^^ = -0.67, pj^ = -0.23, pjg = -0.57. 

The systematic uncertainties, summarized in Table IHl are dominated by uncertainties 
in the backgrounds. Inaccuracies in the simulation can affect the reconstructed neutrino 
momentum, thereby shifting the expected amount of peaking background relative to the 

— K'i^u yield, and hence the extracted —>■ 7r~£^z/ yield. To study such effects, we 
adjust variables in the simulation {Ki production, tracking efficiency, track parameters, and 
shower energy resolution). The sizes of these variations are guided by independent studies 
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FIG. 2: Distributions in (p' for — > K~t^v (left) and —>■ iT~t^v (right), after correcting for 
reconstruction efficiency and smearing in (p. The data include statistical and systematic uncer- 
tainties and are overlaid with various predictions [6-10]. 



TABLE I: The bin yields, with statistical and systematic uncertainties, after correcting for 
efficiency and smearing across bins and normalizing their sum to unity. 

"^(GeV^) ^T/T{K-i+v) ^T/T{-k-1+u) 

[0, 0.75] 0.654 ± 0.010 ± 0.005 0.45 ± 0.05 ± 0.03 

[0.75, 1.5] 0.323 ± 0.015 ± 0.006 0.26 ± 0.06 ± 0.04 

> 1.5 0.024 ± 0.008 ± 0.006 0.29 ± 0.05 ± 0.02 



of the detector and the scale of the small discrepancies observed between data and simulated 
distributions in mhu mhw and Pu- Biases in the simulation can also affect the 7r"£^z/ 
and — > K~£~^i> efficiency ratio and distributions. In practice these effects are small 
since the same selection criteria are applied to both modes and remaining differences depend 
primarily on the decay kinematics, which are readily simulated. We find a small contribution 
from the uncertainties in the efficiencies for successfully identifying hadrons (tt or K) and 
leptons. 

Hadron misidentification, particularly mistaking a kaon from D° K~l'^v for a pion 
from — > TT't^v, poses a serious problem. The probability of misidentifying a kaon as a 
pion is measured as a function of momentum with a sample of K~7i~^ decays. The 

momentum-averaged misidentification probability is (1.9 ± 0.1(stat.))%. We test for differ- 
ences in the misidentification probabilities between kaons from — > K~7r~^ (where a tight 
mass cut is applied) and kaons in our sample (where the mass cut is very loose) by apply- 
ing our technique for measuring misidentification probabilities to simulated events of both 
kinds, and see no hint of bias. However, we see run-to-run variations in the misidentification 
probability that approach statistical significance, and accordingly assign it a conservative 
20% relative systematic uncertainty. 

Additional uncertainty arises from the statistical uncertainty in the — > K~£'^u normal- 
ization, since it determines the background level for 7r~£"'"z/. The branching ratios of 
other semileptonic modes, D° —>■ Xi'^u, relative to D° —>■ K'i'^u also affect the yields, as do 
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the form factors for D K*ii', the charm fragmentation parameters, the background from 
candidates in which a hadron is mistaken for an electron or muon, and the normahzation of 
residual BB events. 

TABLE II: The percent uncertainties in Rq and the normaUzed raw q"^ bin yields. Entries are 
explained in the text. Systematic uncertainties apart from the — > K'£^v normalization and a 
portion of the simulation uncertainty (first row) are correlated between the vr and K modes. 



Unc. Source 






^2" 


^3 


K K K 

af cjf 


Simulation^ 


2.9 


3.4 


4.1 


6.0 


0.5 0.5 2.6 


Part. ID (e)t 


1.9 


0.7 


0.8 


1.3 


0.3 0.3 0.4 


Part. ID (^)t 


2.0 


0.9 


0.9 


1.0 


0.4 0.4 0.6 


K/tt mis-IDt 


3.9 


3.1 


3.5 


3.2 


0.0 0.0 0.0 


Kei' norm.* 


1.0 


1.1 


1.5 


1.9 




Kfiu norm.* 


3.2 


4.8 


3.9 


4.8 




Xli^ br. fr.t 


3.5 


0.8 


1.3 


1.2 


0.3 0.2 0.3 


K*eu f. f.t 


1.1 


0.7 


0.1 


1.7 


0.1 0.1 0.4 


cc frag.^ 


1.7 


0.4 


0.5 


0.0 


0.0 0.0 0.0 


e mis-ID^^ 


0.4 


0.2 


0.5 


0.4 0.0 0.0 0.0 


mis-ID^" 


3.0 


0.9 


0.6 


0.1 


0.0 0.0 0.0 


BB norm.^ 


0.2 


0.2 


0.1 


0.2 


0.0 0.0 0.0 


Total (e) 


6.7 


4.9 


5.8 


7.5 


0.7 0.7 2.7 


Total (fj.) 


8.0 


6.9 


6.9 


8.6 


0.7 0.7 2.7 



Stat, (e)* 7.7 8.0 10.3 15.8 1.3 1.3 3.1 
Stat. (^)* 17.0 24.6 21.9 28.7 4.8 3.7 4.3 

* Assumed uncorrelated across bins, 
f Assumed correlated across q^ bins. 



:|:Assumed correlated across q bins in calculating Rq and uncorrelated across q bins in the form 
factor fits. 

Combining Rq^ and Rq^ after applying a +1% correction to Rq^ to account for the reduced 
muon phase space, gives 

Rq = 0.082 ± 0.006 ± 0.005. 

This result is consistent with the previous world average of 0.101 ± 0.018 j^, but is more 
precise. 

We next determine parameters describing the form factors by fitting the corrected 
distributions. We first use a simple pole parameterization, 

and vary the value of mpoie, constraining the integral over to unity. The quality of the 
fits is good. Dominance by a single pole would imply m^~^^ = ttid*,^. We find m^'^'" = 

l-SGtofatom GeV and m^-^/ = 1.89±0.05lg:[|^ GeV, where the uncertainties are statistical 
and systematic. We also fit the data with a modified pole distribution |lll |. 



{I - / m%. Ml - aq^ / m%,)' 
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to obtain the parameter a. We find a^^"" = 0.37^[}:1? ± 0.15 and a^^^ = 0.36 ± 0.10^[}:S?. 
Our results for m^~^^ and a^^^ suggest the existence of contributions beyond the pure 
pole to the D° —>■ K'i'^u form factor. For — >■ 7i~i'^h', rrip'^'" is consistent with the D* 
mass, though the precision is sufficient to rule out only large additional contributions. 

Several predictions for the form factors are superimposed on our data in Figure |21 Most 
are in satisfactory agreement with the data. The ISGW2 model 0, however, predicts a 
distribution for —>■ K~i~^i' that peaks lower than the data, and accordingly the with 
our data is poor (18 for 2 degrees of freedom). 

Using the value of Rq and parameterizing the form factors with the results of the modified 
pole fit, we find 



+0.006+0.005 

Ktn\\2\-[/ 12 ~ "-'•"-'^'-'-0.007-0.003' 



0.038 



l/f(o)PIK 

where the uncertainties are statistical (±0.003 from Rq and ±0.006 from a) and systematic 
(±0.002 from Rq and to'xm from a). The result is the same within 1% if we use the 
simple pole form factor instead. Using iK^/KsP = 0.052 ±0.001 gives |/^ (0)|/|/:^(0)| = 
0.86 ±0.071 

ao4 i 0.01, where the first error is statistical, the second is systematic, and the 
third is from the CKM matrix elements. This value is consistent with most expectations for 
SU(3) symmetry breaking [7-9,12]. 

We have presented a new measurement of the ratio of n^i^u to —>■ K~£^i' 

decay rates. This result is more precise than the previous world average by nearly a factor 
of two. Our data also provide new information on the K~[^v form factor, a first 

determination of the —>■ form factor and the first model independent constraint 

on |/+(0)||V^rf|/|/:|!^(0)||Vcs|. Together, these offer new checks of SU(3) symmetry breaking 
and the form factors predicted for the semileptonic decays of heavy mesons into light ones. 
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